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ABSTRACT 

Context. In the generally accepted, but poorly documented model, silicate J-type C-stars are binary objects for which the silicate 
emission originates from a circumbinary or a circumcompanion disc. 

Aims. We aim at testing this hypothesis by a thorough spectral and spatial observational study of one object: IRAS 18006-3213. 
Methods. We obtained, analysed and modeled high spatial resolution interferometric VLTI/MIDI observations on multiple baselines 
ranging from 45 m to 100 m. 

Results. All observations resolved the object and show the very compact nature of the N-band emission (~ 30 mas). In addition, the 
highest spatial resolution data show a significant differential phase jump around 8.3 /im. This demonstrates the asymmetric nature of 
the N-band emission. Moreover, the single telescope N-band spectrum shows the signature of highly processed silicate grains. These 
data are used to confirm the model on silicate J-type C-stars for IRAS 18006-3213. We show that the most favourable model of the 
dust geometry is a stable circumbinary disc around the system, seen under an intermediate inclination. 

Conclusions. The data presented on the silicate J-type C-star IRAS 18006-3213 provide evidence that the oxygen rich dust is trapped 
in a circumbinary disc. The formation of this disc is probably linked to the binary nature of the central star. 
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1. Introduction 

The 12 C content in the atmosphere of asymptotic giant branch 
(AGB) stars is believed to increase during the evolution along 
the spectral sequence M — » S — > C. This increase is the result 
of mixing of He-burning material (mainly 12 C and s-process el- 
ements) by the convective envelope through the third dredge-up. 
This mixing process is thought to occur during the relaxation 
after a thermal pulse. The recurrent operation of this dredge-up 
leads to the creation of a carbon star, when carbon becomes more 
abundant in number than oxygen. 

The spectral class ification of carbon stars 
dWaTlerstein & Knappl I1998L and references therein) revealed, 
however, that different spectral groups could be identified and 
only the C(N) stars are now believed to be AGB stars, which 
are in the process of active dredge-up evolution. Other spectral 
classes are much less understood. 

The J-type carbon stars are genuine carbon stars (C>0), but 
they show very low 12 C/ 13 C values, at or near the CNO cy- 
cle equilibrium value of about 4. Moreover, J-type carbon stars 
show n o s-process overabu ndances and often large Li abun- 
dances ( Abia &~Isernl l2000l and references therein). All these 
chemical characteristics are at variance with the AGB chemi- 
cal models and it is fair to say that the class of J-type C-stars are 
chemically far from understood. They do, however, represent a 
significant fraction of all C-st ars (~ 10%) as was shown from a 
recent survey of 1497 stars in lMorgan et alj d2003l) . 

Because of the stability of the CO molecule, O-rich stars cre- 
ate circumstellar envelopes with silicate dust emission, while C- 
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rich stars show C-rich dust features. The detection by IRAS/LRS 
of O-rich dust around C-rich photospheres came therefore as a 
huge surprise dLittle-Mareninll 19861: IWiilems & de Jond [l986). 
Even more intriguing is the observation that all C-stars show- 
ing s ilicate dust are spectroscopically of J-type dLlovd Evansl 
1991). No photospheric distinction can be made within the class 
of J-type C-stars, betw een the silicate emitters and the others 
dOhnaka & Tsui il 19991) . The J-type C-stars showing silicate dust 
emission ar e defined as "silica te J-type C-stars" and constitute 
about 10 % dLlovd Evanslll 99 lh of the J-type C-stars. 

Despite 20 years of research, the evolutionary channel(s) 
leading to the formation of silicat e J-type carbon stars remain(s) 
a mystery. iLittle-Mareninl (|T986) suggested that silicate J-type 
C-stars can be explained by assuming they are binaries con- 
sisting of a C-rich and O-rich giant. However, no direct detec- 
tion of thi s companion has ev er been made, neither spectro- 
scopi callv dLambert et alJll990h nor with speckl e inter ferome- 
try dEngels & Leinertlll994l) . IWiilems & de Jond d 19861) on the 
other hand suggest that these objects experienced a recent tran- 
sition from O-rich to C-rich by means of a thermal pulse. The 
O-rich matter in this hypothesis is then the remnant of a previous 
mass-loss phase, when the giant was still O -rich. This was crit- 
icized from an evolutionary point of view. iLlovd Eva ns (1990) 
argues that such a transition object should only be observed dur- 
ing a few decades, while some of the silicate C-stars are known 
to have C-ri ch photospheres for ove r 50 years. Furthermore, it 
is shown in Yamamura et al.l d2000t) that the silicate feature of 
V778 Cyg did not change in the 14 years time interval between 
the IRAS/LRS and the ISO/SWS observation. Therefore, also 
this scenario seems unlikely. 

The general consensus now about these objects is that they 
are binaries with an unseen companion, probably a main se- 
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quence star dMorrislll987t iLlovd Evanslll990t [YamamuraetaD 
2000). While the star was an O-rich giant, the mass loss was 
captured and stored in a disc, either around the companion or 
around the system. Subsequently, the primary experienced (a) 
thermal pulse(s) and evolved into a C-star. This scenario still 
does not account for the chemical J-type nature of the silicate 
C-star. 

Evidence of this hypothesis has been gathered over the 
years: The long-lived reservoir has been inferred from very 
narrow rotational CO line profiles by Kah ane et al.l d 19981) and 
lJura & Kahand (Il999l) . indicative of Keplerian rather than out- 
flow velocities. Radial velocity measurements of two silicate C- 
stars were consistent with the motion in a bin ary system. No 
concl usive orbit could be detected, however dBarnbaum et al.l 
[1991 . A recent observation of the water maser emission toward 
V778 Cyg at high angular resoluti on suggests the presence of a 
rotating disc dSzczerba et al.l l2~006). From the angular separation 
between this emission and a previous optical observation, the 
authors conclude that this (probably) p rovides support for th e 
circumcompanion model for this object ( Yamamura et al. 2000). 

Despite the more general consensus of the binary nature of 
silicate J-type stars, conclusive direct evidence on their binarity 
is still lacking. It remains also unclear whether one formation 
channel does applies for all silicate J-type C-stars. An open ques- 
tion is also whether the scenario developed for the silicate J-type 
C-stars, could apply for all J-type stars. Whether a possible bi- 
nary nature can have a direct impact on the internal chemical 
evolution of one of the components is yet another open question. 

To investigate the silicate J-type star problems, we studied 
IRAS 18006-3213 at high angular resolution using the MIDI in- 
strument. IRAS 18006-3213 (CGCS 3935, Fuen C 2 70, FJF 270) 
was id entified as a C-star more than 25 years ago by Fuenmavorl 
d!981l) . Soon after, it was realised bv lWillems & de Jongldl986l) 
that this object belongs to the peculiar class of C-rich objects 
showing an O-rich circumstellar environ ment. The J-type na - 
ture of the C-rich giant was confirmed by ILlovd Evans! (11990). 
IRAS 18006-3213 is therefore a perfect candidate to provide an- 
swers to the questions posed above. The distanc e to this object 
is uncertain and was quantified to be 2.6 kpc by Engels (1994) 
based on the matching between the measured energy emitted in 
the infrared, and the assumption that the luminosity is 5000 L , 
which is typical for thermally pulsing AGB stars. 

The paper is organized as follows. In sect. [2] we present the 
observations for which an analysis is given in sect. [3] The data 
provide constraints on the source morphology which are dis- 
cussed in sect. [4] The conclusions are given in sect. [5] 



2. Observations and Reduction 

2.1. observations 

Five interferometric observations of IRAS 1 8006-32 1 3 were per- 
formed using the MIDI instrument at different physical base- 
lines with projected baselines (PB) ranging from 40 m to 102 m. 
One of the observations (with a PB of 40 m) shows bad over- 
lap between both beams and is therefore discarded. A log of 
the successful observations is presented in Table Q] Each ob- 
servation was performed using a prism (with spectral resolution 
A/AA ~ 30) producing spectrally dispersed fringes over the N- 
band (8-13 /im). T he standard observin g sequence was used, 
for which we refer to lLeinert et al.l d2004l) . To correct for instru- 
mental and atmospheric loss of coherence, a calibrator of known 
diameter is measured using the same setup. The visibility cali- 



Table 1. The log of the observations of IRAS 18006-3213. For 
each observation, the corresponding visibility calibrator is given 
as well. 



night 


telescope 


time 


PB 


PA 


visibility 


yyyy/mm/dd 


combination 


hh:mm 


(m) 


n 


calibrator 


2004-04-08 


UT2 - UT3 


08:17 


46.5 


30 


HD 165135 


2004-06-27 


UT1 - UT3 


03:15 


102.4 


25 


HD 134505 


2004-06-27 


UT1 - UT3 


06:06 


96.0 


43 


HD 139127 


2004-06-28 
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08:48 
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49 
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Fig. 1. Observed spectra of IRAS 18006-3213. The MIDI spec- 
trum is plotted using a solid line while the IRAS/LRS spectrum 
is shown using crosses. Clearly the shape in the N-band region is 
very similar between both spectra although the absolute level is 
not identical. Overplotted using a dotted line is the silicate J-type 
C-star V778 Cyg which is shifted by 20 Jy for the sake of clarity. 
The 10 /an profile shows that the silicates around IRAS 18006- 
3213 are much more processed than the dust around V778 Cyg. 



brators of the different science observations are given in Table Q] 
as well. 



2.2. spectroscopic reduction 

A traditional single telescope spectrum is obtained, together with 
each interferometric observation. We corrected for atmospheric 
transmission and obtained flux calibration using the photometric 
observations of the calibrator stars. Theoretical N-band spectra 
of the calibrator stars w ere synthesized using the marcs code 
dGustafsson et alj[l975l and further updates), including the SiO 
opacity. Temperature, gravity and angular diameter of the cali- 
brators were derived from SED fitting by Verhoelst (2005) (cf. 
http://www.ster.kuleuven.ac.be/~tijl/MIDIxalibration/mcc.txt). 
The spectrum of IRAS 18006-3213 obtained for 2004-04-08 
is shown in Fig. Q] Overplotted in this figure is the IRAS low 
resolution spectrum scaled on the basis of the IRAS 12 /mi 
photometry flux. Also shown in this figure is the ISO/SWS 
spectrum of V778 Cyg. 



2.3. interferometric reduction 

To reduce the interferometric data, two different methods were 
employed. The first method is based on power spectrum analy- 
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sis and was executed with the MIA package (MIDI interactive 
analysis, version 1.3) implemented at the Max-Planck Institut 
fur Astronomy in Heidelberg. The second method reduces all 
frames to the same optical path difference (OPD) and was per- 
formed using the EWS package (Expert Work Station) devel- 
oped by W alter Jaffe. For the de tails of both redu ction packages, 
we refer to lLeinert et al.l d2004l) and lJaffd ((2004) respectively. 

In Fig. [2] we show the visibilities obtained with the EWS 
method. The error bars consists both of systematic errors and 
noise type errors. Because the systematic contribution to the er- 
rors is dominant, the actual point-to-point error on the visibility 
is much smaller than indicated in Fig. [2] i.e. the relative visibility 
between e.g. 8 and 9 pm is much better determined than the nu- 
merical value at either wavelength. The visibilities obtained with 
the MIA reduction scheme give very similar results. In the fol- 
lowing, we will work with the EWS reduced data because this 
method has the advantage that the relative phase between the 
different spectral channels is preserved in the reduction. All first 
order information on the phase delay is lost because it is neces- 
sary to compensate for OPD fluctuati ons and firs t-order disper- 
sion effects before co-adding the data (Jaffe 2004). Second order 
effects remain in the data and they can be determined very accu- 
rately. An accuracy of 1 degree RMS phas e error can be ach ieved 
during stable environmental conditions dTubbs etalJ 120041) . In 
Fig. |U we show the instrumental phase of the different observ- 
ing nights as measured for all calibrator objects observed during 
the night. The instrumental phase is clearly very stable with an 
accuracy for all nights of about 5 degree RMS phase error. The 
curvature in the instrumental phase is a result of the varying val- 
ues of the index of refraction of air in the VLTI tunnels. The 
instrumental phase is subtracted from the science target obser- 
vations, largely removing this curvature. As a conservative error 
on the relative phase, we took the RMS phase degree error on 
all calibrator observations during the night. The resulting phase 
information on IRAS 18006-3213 for the different observations 
is shown in Fig. [3] 
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3. Analysis 

3.1. spectral energy distribution 

To construct the spectral energy distribution (SED) of 
IRA S 18006-3213, we used o ptical photometry from the litera- 
ture dLlovd Evansll 1 990i 1 1 99 it) and infrared data from 2 MASS, 
MSX and IRAS. The data are shown in Fig. fusing crosses. The 
SED shows a clear IR-excess representing an amount of energy 
almost equal to that emitted in the optical. 

We note that the DENIS database lists a K magnitude of al- 
most 0.7 + 0.1 mag lower than the one reported in the 2 MASS 
database whilst having compatible J magnitudes. This could in- 
dicate that the (dust) excess emission in the infrared is variable. 
This is, however, incompatible with both the spectral and pho- 
tometric stability of the mid-IR emission (see sect. 13.21 ). In what 
follows we did not take into account the DENIS K band magni- 
tude. 

In Fig.g] the SED of IRAS 18006-3213 is compared to the 
SED of T Lyr, a J-type C-star of similar spectral type but with- 
out silicate emissio n. Both objects are among the reddest m em- 
bers of the class dLlovd Evansl 1 99(jt iLorenz-MartinsIl 996 ) . The 
spectral types indicate low effective temperatures in the range 
of T e ff = 2000 - 2500 K. While both objects show compara- 
ble optical colours, in the IR the SED is completely different. 
IRAS 18006-32 13 shows a large IR excess dominated by thermal 
emission of dust, while the IR spectrum of T Lyr is dominated by 



Fig. 2. The calibrated visibilities for the different measurements 
of IRAS 18006-3213. In the top panel the visibilities for the 
short baseline configurations are shown while in the bottom 
panel those for the longer baseline settings are given. The er- 
ror bars are the RMS value of the different scans as estimated by 
the EWS method. 



molecular features and has a very differe nt spectral slope. The 
total m ass-loss of TLyr is quantified by Bergeat & Chevallier 
d2005l) to be 8.10 _7 M o yr _1 and there is no strong evidence for 
efficient dust production in the wind. 

Since both TLyr and IRAS 18006-3213 have very similar 
spectral type, and have very similar colours in the optical, we 
can assume that also the photosphere of IRAS 18006-3213 does 
not experience a large dust-driven spherical C-rich outflow and 
that the IR spectrum of the underlying carbon star is very sim- 
ilar to the template TLyr. In general, J-ty pe objects show very 
moderate outflows dLorenz-M artins 1993. 

We conclude that the spectral energy distribution of 
IRAS 18006-3213 consists of a photospheric contribution and 
a large infrared excess, likely completely dominated by oxygen 
rich dust alone. 
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Fig. 3. The calibrated phase for the different measurements of 
IRAS 18006-3213. Note that only second order effects in the 
phase variation over the wavelength band can be retrieved. The 
same line styles were employed as the ones used in Fig. [2] Note 
that for the two short baseline measurements (solid line and 
dashed line) the phase is centered around zero while for the long 
baseline measurements (dotted and dashed-dotted line, a large 
phase jump is observed around 8.5 fim. The error bars are the 
RMS values on the instrumental phase during the night. 



Fig. 4. The mean instrumental phase for the different observing 
nights. The curvature of the instrumental phase is caused by the 
varying values of index of refraction of air in the VLTI delay 
line tunnels. The error given is the RMS error on the different 
calibrator observations during the night. The data for 2004-04- 
08 is shifted down by 10 degree while the data for 2004-06-28 
is shifted up by 10 degree for the sake of clarity. Clearly the 
atmospheric conditions were very stable during the observations 
giving rise to very small RMS errors of the order of only a few 
degree. 



3.2. mineralogy 

In Fig. [Q the MIDI spectrum of IRAS 18006-3213 is shown to- 
gether with its IRAS/LRS spectrum. Both show good agreement, 
which indicates that the silicate emission must have been stable 
for over 20 years. This is consistent with the stable emission 
found for the two other silicate J-type C-stars of which the IR 
excess was observed at different epochs (V 778 Cyg and Hen 38, 
Yamamura et al. 2000; Ohnaka et al. 2006, respectively). 

The band profile on the N-band silicate feature is a good 
tracer for the mineralogy of the dust emission. Comparing the 
spectrum of IRAS 18006-3213 with that of V778Cyg (also 
shown in Fig. |TJ, another silicate J-type C-star, it is striking that 
the silicate emission profiles are very different. The profile of 
IRAS 18006-3213 has the signature of highly processed silicates 
(i.e. crystalline and/or large grains) while V778 Cyg shows ev- 
idence for an unprocessed silicate environment (i.e. amorphous 
small grains). 

In order to quantify the difference in mineralogy and sizes 
of the emitting dust particles around both objects, we compute 
a fit through the N-band spectra using calculated emissivities of 
irregularly shaped, chemically homogeneous dust grains. We ig- 
nored the contribution provided by their photospheres because 
this contribution and the effect on the slope of the N-band spec- 
trum is marginal. The most important dust species causing spec- 
tral signature in the 10/mi window are amorphous (i.e. glass 
with olivine stoichiometry) and crystalline olivine, amorphous 
(i.e. glass with pyroxene stoichiometry) and crystalline pyrox- 
ene and amorphous silica. In the fitting, we added also a con- 
tinuum distribution which accounts for the emission by large 
grains and/or for the possible presence of featureless compo- 
nents such as metallic iron, iron sulphide or carbon grains. This 
continuum contribution is modeled using a constant mass ab- 
sorption coefficient. For a reference to the method which is used 



and the assumptions m ade, we refer to iMin et al .] d2005h and 
IVanBoekeletail (120051) . 

In the 10 fim region we are mainly sensitive to the dust grains 
smaller than a few fim. We represented the size distribution of 
the particles by two grain sizes, 0. 1 fim and 1 .5 fim. This size is 
the radius of th e solid sphere w ith an equivalent volume as the 
hollow sphere (Mi n et al ] |2005l) . Particles larger than a few fim 
contribute mainly to the continuum. In addition, we assumed that 
the thermal radiation we analyze originates from optically thin 
regions so that we can simply add the different contributors to 
the spectral emission feature. 

In the fitting process of the N-band spectrum, we assume all 
contributing dust grains, including the ones responsible for the 
continuum, to have the same temperature distribution. Due to 
the limited wavelength range, this temperature distribution can 
be represented by a single Planck curve with one characteristic 
temperature, T c . This is justified because it is likely that the dust 
grains of different species are coagulated, implying thermal con- 
tact between the various components. The characteristic temper- 
atures used in the modeling are given in Table [2] Note that this 
temperature can not be regarded as the temperature of the whole 
dust excess but is only an approximation of the temperature of 
those particles responsible for the N-band emission. 

The abundances of the dust components are determined by 
using a linear least square fitting procedure with constraints on 
the weights to avoid negative values. The temperature of the 
grains and the underlying continuum is varied from to 1500 K 
until a best fit is obtained. The errors on the resulting fit param- 
eters are determined through a Monte Carlo simulation on the 
spectrum. This method ensures that degeneracies between the fit 
parameters show up as large errors in those parameters. 
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Fig. 5. The SED of IRAS 18006-3213 is shown using crosses 
while the SED of a comparable J-type C-star without silicate 
emission, TLyr, is overplotted using diamonds. Shown using 
a solid line is the ISO/SWS spectrum of TLyr and the MIDI 
spectrum of IRAS 18006-3213. While both objects have similar 
colours in the optical, IRAS 18006-3213 shows a huge IR-excess 
in comparison to T Lyr. IRAS 18006-3213 can not be experienc- 
ing a large spherical C-rich outflow. 

Table 2. The composition and grain sizes of the silicate 
dust in the circumstellar environment around IRAS 18006-3213 
(IRAS 18006) and V778Cyg. The percentages for the differ- 
ent particles are listed in mass fraction. The olivine and pyrox- 
ene particles are amorphous while forsterite and enstatite are 
crystalline particles. Two grain sizes were adopted in the fit: 
small particles are 0. 1 fim while big particles are 1 .5 /im large. 
Also given is the total amount of big particles needed in the 
fit and the overall crystallinity fraction. Clearly the dust around 
IRAS 18006-3213 is much more processed as V778 Cyg, both in 
grain size and crystallinity. 





unit 


value 
IRAS 18006 


value 
V778 Cyg 




K 


400 


± 100 


500 ± 100 


small olivine 


% 


0.1 


±0.1 


61.9 ±3.9 


big olivine 


% 


57.4 


± 1.8 


20.6 ± 6.4 


small pyroxene 


% 





±0 


0±0 


big pyroxene 


% 





±0 


0.1 ±0.1 


small forsterite 


% 


3.1 


±0.3 


1.9 ± 0.3 


big forsterite 


% 


15.5 


+ 0.9 


9.3 ± 1.2 


small forsterite 


% 





±0 


0±0 


big enstatite 


% 


15.7 


± 1.1 


6.3 ± 1.4 


small silica 


% 





±0 


0±0 


big silica 


% 


8.1 


±0.3 


0±0 


summary 


big components % 
crystalline components % 


96.8 
34 


±0.5 

±2 


36.2 + 4.4 
17.4 ± 2.4 



The dust parameters and its respective errors for 
IRAS 18006-3213 and V778Cyg are given in Table |2] The 
resulting best fit spectrum is shown as a grey line in Fig. [6] 

The circumstellar silicates around both objects are clearly 
very different. The warm silicates of IRAS 18006-3213 are much 
more processed than those of V778 Cyg. While the fraction of 
large grains is almost 100% (their contribution to the spectrum 
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Fig. 6. The N-band spectra of V778 Cyg (top) and IRAS 18006- 
3213 (bottom) are shown using a black solid line. Overplotted 
using a grey line are the best fit models. Note the almost per- 
fect agreement between the fit and the spectrum. Also shown 
in the figure are: the continuum contribution, the emission pro- 
duced by 1.5 jum sized particles (i.e. the big component), the 
crystalline component and the amorphous component. The dust 
around both objects is clearly very different. While the spectrum 
of V778 Cyg is reproduced mostly through features associated 
with small grains, the circumstellar environment of IRAS 18006- 
3213 consists almost entirely of large (processed) grains. 



is shown using a dotted line in Fig. [6]), they represent only one 
third of the mass for V778 Cyg. Moreover, also the crystallinity 
fraction for IRAS 18006-3213 is twice as large as the one deter- 
mined for V778 Cyg. 

The presence of the high processing observed for 
IRAS 18006-3213, both with respect to grain size and crys- 
tallin ity degree, is indicative of dust stored in a stable disc 
(e.g. lMolster et al.ll2000t IVan Boekel et al.ll2005l) . In O-rich out- 
flows observed around mass-losing objects , the dust is mostly 
in the form of amorphous sili cates (e.g. [Sloan & Price 1995; 
IWaters et alJl 19961: ICamill2002l) . The spectrum of V778Cyg is 
compatible with the emission typically observed in such out- 
flows. We can conclude that the degree of processing of the cir- 
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Fig. 7. The angular size of the emission in the N-band as deter- 
mined using a uniform disc model. The fit was performed using 
the two visibility measurements at the shortest baselines only 
and should therefore only be considered as a rough first order 
estimate. Note that the point-to-point error is much smaller than 
the error given in the plot. 



cumstellar environment of silicate J-type C-stars is clearly not 
uniform. 



3.3. N-band visibilities 

In this section we discuss the global characteristics of the N- 
band emission as constrained by the visibilities shown in Fig. [2] 
The object is clearly resolved at all baselines. In addition, a 
strong wavelength dependence is observed in the visibility ob- 
servations. 

To quantify the global characteristics of the N-band emis- 
sion, we modeled the visibility measurements assuming using a 
simple geometry, i.e. a uniform disc intensity distribution. We 
limit this part of the analysis to the two visibility curves at the 
lowest spatial frequency, because only these probe the large scale 
structure. The highest spatial frequency data yield very low vis- 
ibilities and their interpretation is therefore very sensitive to the 
exact unresolved photospheric contribution. The possible pres- 
ence of small scale structures in the dust excess like clumps, 
warped structure of the disc, etc.. are very poorly constrained 
with the limited number of visibilities observed. The long base- 
line measurements are, however, used for their phase informa- 
tion, see sect. 13.41 

The resulting angular diameter of the N-band emission is 
shown in Fig. Q The dust emission is clearly compact (25 mas 
at 8 fim to 35 mas at 13 yum) and it shows a clear wavelength de- 
pendence. Globally, the size increases as wavelength increases. 
This is likely the result of a temperature distribution, with the 
colder dust located in a larger structure than the hotter dust. In 
addition, a clear dip in angular size is observed starting around 
9 fim. Since the error bars in Fig.Qrepresent mostly a scaling un- 
certainty in the fluxes, the relative difference observed between 
e.g. 8.5 fim and 9.5 fim is significant. This dip in the derived an- 
gular size is the result of the shape of the visibility curve, which 
increases rather steeply around 9 fim and flattens out till 11.5 
fim after which the visibility decreases (see top panel in Fig . [2j) . 



Note that a Gaussian intensity profile provides the same qualita- 
tive results. 

Given the amount of emission and the size of the emission 
region at 8 fim, which originates from continuum emitting par- 
ticles (see Fig. it is clear that the continuum emission must 
come from an optically thick region, i.e. r > 1: If the emission 
were to come from an optically thin region, the derived tem- 
perature would be inconsistent with the SED. For example, if 
t = 0.1, then a temperature > 1200K is necessary to reproduce 
the flux at 8 fim. 

IR AS 18006-3213 was identified as a C-star by iFuenmavorl 
(1981). Therefore, any O-rich dust present in the circumstellar 
environment must have been ejected by the AGB giant at least 
25 years ago. Assuming a typical outflow velocit y of lOkms" 1 
and a distance for IRAS 18006-3213 of 2.6 kpc ( fengelslll994l ). 
the O-rich dust would need ~ 15 years to arrive at the derived 
radius of the dust emission. This shows that the current size of 
the disc is in contradiction with an outflow scenario, since the 
actual central carbon star will not produce silicate dust. 

The fact that the "feature" in the visibility profile has the 
same shape as the silicate emission profile (see Fig. [TJ leads us 
to believe that the dip in the angular size comes from a silicate 
emission feature originating fr om a smaller r adius than the un- 
derlying continuum emission. Ohnak aet alj d2006l) observes a 
similar trend for the silicate carbon star IRAS 08002-3803 and 
they show from a careful modeling that a two-species model can 
reproduce the observations. In their model, the silicate compo- 
nent has a steeper temperature gradient than the other component 
due to the different opacities of both particles. Therefore, the 
emission from the silicate component originates from a smaller 
angular size. However, this is not the only possible interpreta- 
tion. An alternative explanation for the observed "feature" in the 
visibility curve is that different grain types do have a different 
spatial density distributions. This was shown alrea dy in some 
discs around Herbig Ae stars (van Boekel et al.l2004 ) and for the 
disc a round the binary post- AGB object HD 52961 ( Deroo et all 
l2006h . 

In summary, the N-band interferometric observations clearly 
resolved the circumstellar environment of IRAS 1 8006-32 1 3 and 
revealed the compact nature of the emission region. The size of 
the emission provides strong indication that the N-band emis- 
sion must come from a compact stable geometry in the system. 
Moreover, the observations show that the particles responsible 
for the silicate emission feature must be emitting from a smaller 
region than the particles responsible for the continuum emission. 

3.4. N-band phases 

The calibrated phases (Fig. show a large phase-jump around 
8.3 fim for the two observations at the longest baseline. This 
phase jump does not equal 180 deg and/or it is not accompanied 
with a visibility transition through zeroj. Moreover, the smooth 
shape of the phase curvature indicates that it is a slowly evolv- 
ing feature. This phase information shows that the source is not 
centrally symmetric with an identical photocenter throughout the 
N-band. 

The photospheric contribution at 10 fim of the total received 
flux is 5 ± 1 % (assuming a blackbody with the expected T e g , see 
sect. 13. U . At the longest baselines, the dust emission is strongly 



1 Note that, due to the finite bandpass of MIDI, measuring a visibil- 
ity which equals zero is impossible. In reality, a monochromatic zero 
visibility would correspond in our polychromatic measurement with 
V~ 0.008, which is clearly not observed in the measurements. 
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resolved. The unresolved photosphere therefore becomes a sig- 
nificant contributor to the correlated flux of the system. The 
asymmetry in the system is therefore most likely explained as- 
suming an offset between the photocenter of the dust excess and 
the giant photosphere. 

A photocenter offset between both components (the unre- 
solved photosphere and the resolved circumstellar dust emis- 
sion) results in the observed phase through the combination of 
the following effects: 

- The components, with different positions on the sky, have 
different spectral energy distributions. Therefore, the flux of 
both sources will contribute different fractions of the total re- 
ceived flux at different wavelengths, producing a phase dif- 
ference between the spectral channels. 

- The components have different positions on the sky and this 
displacement is more resolved at 8 fim than at 13 //m, leading 
to a modulation of the phase signal with the spatial resolu- 
tion. 

The probability that the phase-jump is explained solely by 
the second effect is very small, since the same phase-jump is 
observed at different instrumental setups (see Fig. [3j. Indeed, 
finding the same shape in phase with different projected base- 
line lengths and angles can only be understood if we assume a 
very specific source geometry: the projected angle and projected 
asymmetry of the source on the sky must be such that they cancel 
in the different observed setups. The probability of this geometry 
is clearly extremely low. 

To see how the shape of the phase over the N-band can be re- 
produced by different spectral signatures of both components (in 
combination with the changing spatial resolution over the wave- 
length band), we quantified a simple illustrative model. 

The most important non-linear feature in the phase is cen- 
tered around 8.3/im (see Fig. (3). For the observations at the 
longest baselines, a strong positive slope is observed before 
8.3 fim. Between 8.3 fim and 9.3 fim the phase strongly declines. 
In our model, an absorption/emission feature in the giant photo- 
sphere and/or the dust spectrum could cause these slopes. The 
most straightforward identification is the silicate feature in the 
dust emission itself. This feature changes the flux ratio drasti- 
cally between 8.3jum and 9.3jum (see Fig.[6]l. While the emis- 
sion from the dust excess increases by almost 50 % over this 
spectral region, the emission of the giant photosphere decreases 
by 20% in the assumption of a Rayleigh- Jeans law. Clearly, 
the global photocenter will shift rapidly toward the photocen- 
ter of the dust emission moving from 8.3 fim to 9.3 fim. In the 
spectral region before 8.3/vm, the presence of a strong absorp- 
tion/emission feature in the dust spectrum is excluded. The un- 
resolved photosphere is the only component for which a large 
change in flux contribution over that specific spectral region 
would go unnoticed in the global N-band spectrum. The most 
likely source is the presence of a molecular photospheric absorp- 
tion band due to C2H2 and HCN located around 7.5 //m. This fea- 
ture appears to be omnipresent in C-stars as is demonstrated by 
ISO dAoki et al.|[l999t lJ0rgensen et alj|2000t lYang et alJ [2004: 
Gautschv-Loidlet al.ll2004l) and Spitz er dMatsuura et afl [2006 ; 
Sloan et all l2006t IZiilstra et all 120061) . As a template of what 
this molecular band for IRAS 18006-3213 could look like, we 
took the ISO/S WS spectrum of T Lyr, a J-type C-star of similar 
spectral type but without silicate dust. Its photospheric molecu- 
lar band at 7.5 fim is likely somewhat polluted by its moderate 
mass loss, but the band shape is expected to remain the same. 
Moreover, the double-peaked shape of the band of TLyr (see 
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Fig. 8. A model for the correlated flux ratio between the infrared 
excess and the photospheric emission of IRAS 18006-3213. For 
the infrared excess we assumed a visibility of 5 % over the wave- 
length band and for the photospheric emission we assumed the 
photospheric emission of T Lyr as reference. To be able to com- 
pare with the observed phase behaviour over the wavelength 
band (see Fig. 0), we removed the first order behaviour of this 
ratio. 

Fig- EJ is representative for the typical band shap e due to C2H2 
and H CN: see e.g. the left panel of Fig. 1 1 in lLagadec et al.l 
(2006), where the averaged spectra around the 7.5yum band of 
14 SMC C-stars is shown. 

In Fig. [8] we show the correlated flux ratio between both 
components assuming an unresolved photosphere and a constant 
visibility over the wavelength band for the dust excess. In or- 
der to compare directly the wavelength dependence of this ratio 
with the observed phase (see Fig. [3j, we removed its first or- 
der behaviour over the wavelength band in Fig. [8] The observed 
phase shape over the N-band is reproduced well with the phase 
jump located at the observed wavelength position. The strength 
of the phase jump around 8.3 fim is enhanced if we assume the 
infrared dust excess to have a lower visibility. The shape of the 
phase jump observed at the longest baselines is therefore ex- 
plained very well by: A shift of the global photocenter toward 
the giant photosphere in the region between 7.8 fim and 8.3 fim, 
followed by a shift toward the photocenter of the dust excess 
between 8.3 //m and 9.3 fan. 

The low visibilities at the longest baselines are determined 
by a complex addition of an unresolved photosphere and a cir- 
cumstellar dust component with a different phase. This dust 
component has an unconstrained (inner) geometry. To provide 
nevertheless a quantified estimate of the offset between both 
components, we assumed in our model the following: (i) The 
emission of the C-giant can be represented by the photospheric 
emission of T Lyr. ( ii) We assume that the wavelength depen- 
dence of the diameter of the dust excess obtained at the shorter 
baselines (see Fig. 13 is also applicable at the longest baselines. 

Given the previous considerations and assumptions, we mod- 
eled the visibility and phase observed at the 96.0 m baseline us- 
ing two free parameters. The first parameter is a shift down in 
the size of the dust excess as determined from the measurements 
at the shortest baselines (see Fig- [2J, to simulate geometrical ef- 
fects at high spatial resolution. The second parameter is the off- 
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Fig. 9. The inverse^- 2 distribution of the minimalisation between 
our model and the observations at the 96.0 m baseline. Four con- 
tours are shown ranging from 0.15 to 0.30 in steps of 0.05. Two 
solutions are easily identified with distances between the giant 
photosphere and the dust excess of 2.4 mas and 23.5 mas. The 
optimum solution is compared to the observations in Fig.fTUl 



set between the giant photosphere and the photocenter of the 
dust excess. In Fig.|9j the inverse^ 2 distribution is shown. Two 
solutions are found: The optimum solution provides a separation 
between both components of only 2.4 mas (x 1 = 3.1), while the 
other provides a much larger separation of 23.5 mas (x 2 = 3.5). 
The resulting fit to the visibility and phase for the optimum so- 
lution is shown in Fig. [10] Both observables are reproduced very 
well considering the simplicity of the model. Modeling the mea- 
surement at the 102.4 m baseline provides similar results with 
best fit separations between both components of 4.7 mas and 
14.9 mas. We remark that under the assumptions we made, the 
visibilities confirm the non-zero phase which is measured. To 
model the shape and the absolute value of the visibility ampli- 
tude, a complex (i.e. non-real) addition of the visibilities is nec- 
essary. 

4. Discussion 

4.1. location of the dust 

Silicate C-stars are identified as objects displaying a C-rich pho- 
tosphere with an O-rich circumstellar environment. Within all 
the known silicate C-stars, large differences in silicate mineral- 
ogy are observed. This was demonstrated in sect. I3.2I and ap- 
parent fro m the wide variety in 10 /mi silicate profiles (see e.g. 
Fig. 4 in iLlovd Evansl [T990) for silicate J-type C-stars. While 
some objects show a mineralogy which resembles the miner- 
alogy observed in O-rich AGB outflows (e.g. V778 Cyg and 
BM Gem), other show signatures of highly processed grains (e.g. 
IRAS 08577-6035 and EU And). The object which is at the fo- 
cus of this study, IRAS 18006-3213, clearly belongs to the latter 
class and shows a crystallinity degree of 35 % and a fraction of 
large micron sized parti c les of almost 100 %. 

In lYamamura et alj d2000h . it is shown that the geometry 
which best reproduces the ISO/SWS observation of V778 Cyg, 
consists of an O-rich reservoir located around the companion 
star. Because O-rich AGB stars typically show amorphous out- 
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Fig. 10. The observed visibilities and phases at the 96.0 m base- 
line are compared to the best fit model with a separation of 
2.4 mas between the giant photosphere and the photocenter of 
the dust excess. The model reproduces the characteristics of the 
data very well. Note that the worst part of the fit is located before 
8.3 micron, where the effect of the uncertainty in the shape and 
strength of the molecular band in the giant photosphere is large. 



flows, the material captured by the component is likely domi- 
nated by amorphous grains. Once in the disc, part of the material 
in the inner region close to the low luminosity companion, may 
be heated above the glass temperature at which point it will be- 
come crystalline. The processed grains will likely not provide a 
strong signature in the integrated spectrum and the spectrum will 
be dominated by the evaporation of the amorphous outer layers 
of this circum-companion disc. This evaporation is triggered by 
the strong radiation pressure of the primary AGB giant. 

In this scenario it is, however, difficult to account for the high 
processing degree observed for the circumstellar environment of 
IRAS 18006-3213. The alternative geometry is when the circum- 
stellar material is captured in a circumbinary disc. In this geom- 
etry, the inner material facing the central binary, is heated above 
the glass temperature and can become crystalline. Because of 
the high-luminosity AGB giant, this region is significantly larger. 
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The impact of this material on the spectrum is therefore expected 
to be significant, thus providing an overall processed feature. 

The visibility measurements of IRAS 18006-3213 also 
favour a geometry where the dust is stored in a circumbinary 
disc. If the silicate emission were to come from evaporation of 
the amorphous outer layers of the disc, this emission would be 
larger than the underlying continuum. Instead, the opposite is 
observed for IRAS 18006-3213. 

4.2. asymmetry in the system 

As shown in Sect. 13.41 asymmetry is present in the source geom- 
etry. This asymmetry is interpreted as originating from a photo- 
center offset between the contribution of the unresolved C-rich 
AGB photosphere and the resolved thermal emission of circum- 
stellar dust. An offset between both components can be a natural 
result of either: 

- the motion of the AGB giant in its (suspected) binary orbit 

- the presence of a circumbinary disc seen under a non pole-on 
inclination 

- the presence of a disc around the companion 

While these scenarios provide different physical offsets be- 
tween the giant and the dust excess, the expected asymmetries 
are compatible with the solutions found in sect. 13.41 The MIDI 
spectrum and the visibility amplitudes however point to the oc- 
currence of a circumbinary disc rather than a circumcompanion 
disc. We therefore evaluate one of the first two options more 
likely. 

We conclude that - although the interferometric phase infor- 
mation is not conclusive - an inclined circumbinary disc is the 
most probable source of asymmetry in the system, possibly in 
combination with the binary motion of the AGB giant. 

4.3. comparison with binary post-AGB objects 

IRAS 18006-3213 shows remarkable similarities to the binary 
post-AGB objects which are known to be surrounded by a cir- 
cumbinary disc, since the detected orbits are all w ithin the dust 
sublim ation radius (for a recent review we refer to lVan Winckell 
2003). These objects have similar high luminosities and it is 
shown that high processing degrees are observed for a very high 
fraction of objects (over 60 % of t he N-band spect ra show a crys- 
tallinity fraction larger than 40 %. lDe Ruvterl2005l) . Also the vis- 
ibility measurements of IRAS 1 8006-32 1 3 provide a very similar 
size and spectral behaviour as observed for the only b inary post- 
AGB object resolved with MIDI (i.e. HD 52961. fDeroo et alj 
2006). This is shown in Fig. QT| where the absolute angular di- 
ameter is given for both objects. The similarity between both 
curves indicates that similar physical processes are active in both 
objects, providing again evidence for the circumbinary nature of 
the disc around IRAS 18006-3213. We do not expect to see a 
similar phase jump for HD 52961 as for IRAS 18006-3213 be- 
cause the post-AGB photosphere of HD 52961 is much hotter 
(and O-rich), making the photospheric spectrum devoid of strong 
molecular lines. 

5. Conclusions 

The spectroscopic and interferometric data presented on 
IRAS 18006-3213 provide clear and direct evidence that the 
oxygen rich reservoir is located in a compact region close to the 
star: 
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Fig. 11. The angular diameter of IRAS 18006-3213 (solid line) 
is compared to the angular diameter of the dust emission around 
HD 52961 (squares). The latter object is a binary post-AGB ob- 
ject surrou nded by a circumb inary disc with a binary period of 
1300 days dDeroo et alj |2006). Clearly both curves provide a re- 
markable similarity both in size and in spectral behaviour. 

- The stability of the N-band spectrum and the observation 
that the circumstellar environment consists of highly pro- 
cessed silicate grains both indicate the existence of a long- 
lived reservoir in the system where the processing occurred. 

- The interferometric observations show the very compact na- 
ture of the N-band emission, again providing strong evidence 
that this emission originates from a long-lived Keplerian 
disc. 

- The interferometric observations show an asymmetric N- 
band emission. This provides direct evidence of the non- 
spherical nature of the emission. 

Using the spectroscopic and interferometric data, we show 
that the most favourable geometry supported by the data is a 
long-lived disc-like reservoir around the system which is seen 
under an inclination. The formation of this disc is likely linked 
to the strongly suspected binary nature of this object. We ar- 
gued that the bimodal appearance of the silicate N-band profile 
in silicate J-type stars (processed versus non-processed grains), 
could be the result of a different location of the dust in the sys- 
tems which is either circumcompanion (e.g. in V778 Cyg) or cir- 
cumbinary (e.g. in IRAS 18006-32 13). The testing of this will re- 
quire direct detection of the orbits, which will need extensive ra- 
dial velocity monitoring campaigns and/or interferometric mea- 
surements providing (closure) phase information. 
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P. Deroo et al.: The circumbinary disc around the J-type C-star IRAS 18006-3213 
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